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ABSTRACT 

We present an in-depth case study of three molecular clouds associated with the 
walls of the Galactic supershells GSH 287+04-17 and GSH 277+00+36. These clouds 
have been identified in previous work as examples in which molecular gas is either being 
formed or destroyed due to the influence of the shells. 12 CO(J=1-0), 13 CO(J=1-0) and 
C 18 O(J=l-0) mapping observations with the Mopra telescope provide detailed infor- 
mation on the distribution and properties of the molecular gas, enabling an improved 
discussion of its relationship to the wider environment in which it resides. We find that 
massive star formation is occurring in molecular gas likely formed in-situ in the shell 
wall, at a Galactic altitude of ~ 200 pc. This second-generation star formation activity 
is dominating its local environment; driving the expansion of a small Hll region which 
is blistering out of the atomic shell wall. We also find new morphological evidence of 
disruption in two smaller entrained molecular clouds thought to pre-date the shells. 
We suggest that at the present post-interaction epoch, the lifetime of this surviving 
molecular material is no longer strongly determined by the shells themselves. 

Subject headings: ISM: bubbles, ISM: clouds, ISM: evolution, ISM: molecules, radio 
lines: ISM, stars: formation 

1. Introduction 

The last decade has seen the emergence of a new paradigm for molecular cloud formation, 
in which the production of star-forming molecular gas is integrated into the modern picture of 
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a dynamic, turbulent interstellar medium (ISM). In this picture, the compression, cooling and 
fragmentation of the diffuse atomic medium in colliding turbulent flows produces clumpy sheets 
and filaments of cold material, which may go on to become molecular and self-gravitating if density 



and column density conditions are met 


(Hennebelle & Perault 


1999; 


Koyama & Inutsuka 


2000 


Hartmann et al.|2001;|Bergin et al.||2004 


Audit & Hennebelle||2005; Vazquez-Semadeni et al.|2006 


2007 


Heitsch & Hartmann||2008 


Heitsch et al. 


2008, 


Inoue &; Inutsuka 


2009). 



Within the scope of physical systems covered by this theory are supershells - large loops, 
shells and cavities, often hundreds of parsecs in diameter, that are formed by the cumulative stellar 



feedback from OB clusters (e.g. |Heiles||1979} |McCray Kafatos| [1987] |McClure-Griffiths et al. 
2002). Multiple stellar winds and supernovae drive a shock front into the surrounding medium. 



which persists over large enough size-scales and timescales that the transition from atomic to 
molecular gas may occur in the cooled material of the swept-up shell. 



Previous work by Dawson et al. (2011, hereafter Paper I) has presented parsec-resolution 
observations of the atomic and molecular ISM in two Galactic supershells, in which we find an 
enhanced level of molecularization over the volumes of both objects. This provides the first direct 
observational evidence of increased molecular cloud production due to the influence of supershells, 
demonstrating that in the case of these two objects the global influence on the molecular ISM 
is a positive one. However, our observations also highlight another long-recognized fact - that 
the influence of shells may also be disruptive to molecular gas on local scales; driving the atomic- 
molecular transition in the opposite direction in pre-existing clouds disrupted by the shell's passage. 

Here we zoom in on three objects from this previous work in order to present a detailed case 
study of molecular clouds in supershell walls. We carry out observations in the 12 CO(J=1-0), 
13 CO(J=1-0), and C 18 O(J=l-0) lines with the Mopra telescope, and combine these with Hi data 



from Paper I and McClure-Griffiths et al. (2003), to examine objects in which the molecular-atomic 
transition is likely proceeding in opposite directions. The Mopra data reveal the molecular gas at a 
resolution of ~ 45", providing detailed information on its distribution and properties, and enabling 
a greatly improved discussion of its relationship to the wider environment in which it resides. 

For the remainder of this introduction we describe the target objects, giving background de- 
tails on the properties of the molecular clouds and their parent shells. The CO observing and data 
reduction strategy are described in Section [2] Section [3] deals in turn with the three clouds; pre- 
senting molecular gas distributions, deriving masses and column densities, assessing star formation 
activity, and presenting an in-depth examination of each in terms of its relationship to its local 
environment and parent shell. Our findings are summarized in Section |4| 



1.1. Target Objects 



We observe three known molecular clouds in the Galactic supershells GSH 287+04-17 and 
GSH 277+00+36. Both shells are large, evolved objects, with radii of > 150 pc, swept-up masses of 
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> 1O 6 M , ages of > 10 7 yr, and gentle expansion velocities of 10-20 km s 1 (for details see Fukui 



et al. 1999 McClure-Griffiths et al. 2000 2003 Dawson et al. 2008b 2011). Their distances are 



estimated as 2.6 ±0.4 kpc and 6.5 ±0.9 kpc respectively, with GSH 287+04-17 located in the inner 
and GSH 277±00±36 in the outer Galaxy. Crucially, both contain large quantities of associated 
molecular gas, equivalent to ~ 10-20% of their total neutral gas mass. The target regions are 
shown in Figures [T] and [2] An additional fourth region covers a small section of shell wall in which 
no CO has been previously detected. 

Cloud A comprises the lower half of a large (M ~ 7700 Mq) molecular cloud embedded in the 
wall of GSH 287±04-17. At a distance of z ~ 200 pc above the Galactic midplane, its altitude is 
several times greater than the local molecular gas scale height (<r 2 ~ 60 pc, e.g. Malhotra 1994). 



This object was first noted by Fukui et al. ( 1999 ) and later appears in the cloud catalogue of Dawson 



et al. (2008a) as G285.90±4.53 (cloud 16), detected in both 12 CO(J=1-0) and 13 CO(J=1-0). Paper 



|I] identifies this cloud as a likely example of molecular material formed in-situ in the shell walls, 
based on the configuration of the Hi and CO emission, the consistency of formation timescales, 
and the presence of sufficient material for shielding of the CO molecule. The molecular material 
forms an integral and coherent part of the shell wall, with CO emitting gas distributed along the 
curved rim of the atomic shell, and a tight correlation between the Hi and CO emission observed 
throughout the feature. A bright, young star forming region dominates the lower regions of the 
cloud. 

Cloud B is also located in the wall of GSH 287±04-17, and is identified as G288.27±5.60 (cloud 



92) in the catalogue of Dawson et al. (2008a). The shell walls show tapered fingers of Hi pointing 



radially inwards towards the shell interior, some of which harbor clumps of molecular gas at their 
tips. With an estimated mass of M ~ 900 Mq, Cloud B is the largest of these molecular clumps, 
and is located at the very tip of an extended tail of Hi. Unlike Cloud A, this object is not associated 
with large Hi column densities or well-embedded within a section of wall. It is instead offset towards 
the very edge of the Hi finger, with CO emission extending outside the detected limits of the atomic 



feature. Paper I identifies this cloud as a likely example of remnant molecular material that existed 



in the region before the shell's passage, which has subsequently become disrupted and entrained in 
the shell interior. Morphological signatures and a lack of material for shielding of the CO molecule 
attest to this interpretation. This cloud is also a high altitude object; located at a distance of 
z ~ 250 pc from the Galactic midplane. 

Cloud C is a small molecular cloud located at the tip of an Hi 'drip' in the wall of GSH 
277±00±36. This cloud was detected at a single position in 12 CO(J=1-0) in Paper I and has 



not yet been mapped. It's location at the tip of an Hi feature marks it as another candidate for 
pre-existing molecular material that is interacting with the expanding shell. Like Clouds A and B, 
it is also located some distance from the Galactic midplane, at z ~ 170 pc. 

The region labelled 'D' in Figure [2] covers another Hi 'drip' in the walls of GSH 277±00±36, 
in which no molecular emission has previously been detected. Its appearance and location are both 
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similar to the feature that hosts Cloud C. 



2. Observations 

The four target regions were observed with the Mopra telescope, near Coonabarabran, Aus- 
tralia. Observations of the GSH 277+00+36 regions were made in June 2005 and were targeted 
towards positions of high Hi intensity, since no molecular detections had previously been published 
at that time. Observations of the GSH 287+04-17 clouds were made in October 2006 and October 
2009, and covered all positions at which 12 CO(J=1-0) emission had been previously detected with 



the NANTEN telescope (Fukui et al. 1999). 



The transitions observed were the 12 CO(J=1-0), 13 CO(J=1-0) and C 18 O(J=l~0) lines, at 
115.3, 110.2 and 109.8 GHz. In 2005 the Mopra telescope was equipped with a dual band SIS 
receiver that provided 600 MHz instantaneous bandwidth between 86 and 115 GHz. For 12 CO 
observations the correlator was configured to provide 1024 channels over a bandwidth of 64 MHz, 
providing a velocity resolution of 0.16 km s _1 and a velocity coverage ~ 165 km s _1 . For 13 CO 
and C 18 the spectral resolution was doubled to 0.08 km s -1 , reducing the total velocity coverage 
to ~ 80 km s" 1 . The 2006 and 2009 sessions made use of the newly installed MMIC receiver and 
the MOPS digital filter bank, which can simultaneously record dual polarization data for up to 
sixteen 137.5 MHz zoom bands positioned within an 8 GHz window. Zoom bands centered on the 
rest frequencies of the three lines each contained 4096 channels, providing a velocity resolution and 
coverage of 0.09 km s _1 and ~ 360 km s _1 in all lines. The remaining zoom bands were positioned 
to cover other dense gas tracers, but no emission was detected. 

Observations were made in an on-the-fly (OTF) raster mapping mode, in which the telescope 
records data continuously while scanning across the sky. Mosaics of 5' x 5' maps were arranged to 
cover the target regions, with an overlap of 12" between each map to ensure smooth coverage. Each 
map was observed at least twice in orthogonal scanning directions to minimize scanning artifacts. 
For the 2005 session the scan speed was 3 arcsec s _1 and spectral data was recorded at 6" intervals 
along the scan direction. The spacing between scan rows was 8" in 12 CO and 9" in 13 CO and C 18 0. 
In 2006 and 2009 the scan speed was 3.5 arcsec s _1 , the sampling interval was 14" and the spacing 
between scan rows was 10". These settings fulfill the minimum requirements for oversampling of 
the ~ 33" (FWHM) Mopra beam. For all sessions an off-source position was observed once per 
scan row. 

The pointing solution of the telescope was verified once every 60-100 minutes via observations 
of the SiO masers X_Cen, VY_CMa, RWJVel and R_Car, and corrections were applied for pointing 
errors of greater than 10" in either azimuth or elevation. System temperature measurements were 
made with a paddle at the start of each map and at 30 minute intervals thereafter, and took typical 
values of 450-650 K at 115 GHz and 250-350K at 109 and 110 GHz. 



Bandpass calibration, baseline subtraction and calibration onto a T\ scale are performed with 



- 5 - 



the livedata package. The spectra are then gridded into cubes using gridzilla. Q The data are 
weighted by the inverse of the system temperature, and convolved with a truncated Gaussian 
smoothing kernel with a FWHM of 60" and cutoff radius of 30" to improve signal-to-noise ratio. 
This results in a final effective angular resolution of 45". Post-gridding, a 27" band was trimmed 
from the outer edge of each image to remove noisy map edges. 

The gridded data are then converted to a main beam temperature scale using a scaling factor 
determined by daily observations of the standard calibrator source Orion KL (aeigso = 5 h 32 m 47.5 s , 
,5 B1950 = -5°24 / 21"). |Ladd et al.| ( [2005[ ) have shown that for 12 CO(J=1-0) and 13 CO(J=1-0) Orion 
KL peak main beam temperatures recorded with Mopra match those measured by the 15m SEST 
telescope, demonstrating the source is sufficiently extended within the beams of both instruments 
that measurements from the two may be directly compared. Referring to SEST documentation, 
and assuming a SEST main beam efficiency of r/ m b = 0.7 ( Johansson et al.||1998 ), scaling factors of 
rj = 0.53 ±0.02, 0.40 ±0.02 and 0.53 ±0.02 are required to rescale the Mopra 12 CO(J=1-0) data for 
2005, 2006 and 2009, respectively. The first and last of these figures are in excellent agreement with 
the documented Mopra 'extended main beam efficiency', 7? x b ~ 0.55, appropriate for observations 
of sources with angular extent > 80" ( Ladd et al.|2005 ) . The outlying 2006 figure is consistent with 
data recorded by other Mopra observers in the same period ( Hughes et al.|[2010 ). The equivalent 
scaling factors derived for 13 CO(J=1-0) are 0.59 ± 0.02 (2005), 0.59 ± 0.02 (2006) and 0.58 ± 0.03 
(2009). These are also applied to the C 18 O(J=l-0) line. 

In cases where a region is covered by data from two different years, the spectra are re- gridded 
to a common grid and weighted by the inverse of their rms noise before being linearly combined. 
Finally, narrow channel data are binned by two velocity channels, bringing the final velocity res- 
olution of all datasets to 0.16-0.18 km s _1 . The rms noise per channel varies across the observed 
regions, and takes values of 300-600 mK in 12 CO(J=1-0), and 100-200 mK in 13 CO(J=1-0) and 
C 18 O(J=l-0). 



3. Results & Discussion 

12 CO(J=1-0) and 13 CO(J=1-0) emission was detected in Clouds A, B and C, and C 18 0(J=1- 
0) was detected in Clouds A and B only. No CO was detected in Region D, to a 3c detection limit 
of 1.0 K in 12 CO(J=1-0). 



1 Binaries and source code for livedata and gridzilla are available from 
http: / /www. atnf.csiro.au/computing/software/livedata. html 
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3.1. Derivation of Column Densities & Masses 



Before examining each of the three clouds in depth, we briefly outline the common analysis 
performed for each. 

Figure [3] shows H2 column density maps for the three clouds. For material traced in 12 CO(J= 



1- 

0) we use an X-factor of X = 1.6 x 10 2U K -1 km s _1 cm~ 2 ([Hunter et al.||l997|) to convert CO 



integrated intensity directly to molecular hydrogen column density, Nn 2 • This treatment implicitly 
assumes that the 12 CO line is optically thick and that the molecular clouds are self-gravitating (see 
e.g. |Heyer et~aL||200l|>. For 13 CO(J=1-0) and C 18 O(J=l-0), iYi3 CO and N c i So may be obtained 



from the following relations under the assumption of local thermodynamic equilibrium (LTE) : 



iVi3 CO /cm~ 2 = 2.42 x 10 14 
iV c i8o/cm" 2 = 2.54 x 10 14 



T ex Av 



1 - exp(-5.29/T e: 

T ex Av 
1 - exp(-5.27/T e: 



^ Tl3 CO,i ! 

i 



(1) 

(2) 



where Av is the channel width in km s 1 and the sum is performed over all channels, i, within the 
cloud velocity range. Here, the optical depth in the ith channel, Tj, is computed as 



Tl3 co, 



T c 18 o,. 



In 



In 



T m b .13 



CO, 



5.29 ( {exp(5.29/T ex ) - l}" 1 - 0.164 



T m b,i,C ls O,i 



5.27 ( {exp(5.27/T ex ) - l}' 1 - 0.166 



(3) 



(4) 



CO column densities are then converted to Nn 2 by assuming H2-CO abundance ratios of 5 x 10 5 
for 13 CO ( |Dickman|[T978| and 6 x 10 6 for C 18 ( |Frerking et"al|l982 ). 



In the above expressions, the excitation temperature, T ex , is estimated from the peak temper- 
ature of the 12 CO(J=1-0) line, which is assumed to be optically thick at all positions except diffuse 
cloud envelopes. We take T ex ~ 12.5 K for all clouds, based on typical peak 12 CO(J=1-0) main 
beam temperatures in clouds A and B outside localized regions of active star formation. Cloud C 
T m b values are not included in this estimation, since its large distance and small angular extent 
mean that the filling factor of emission within the beam is likely to be low. 

The molecular mass traced in each line is obtained from the sum of iV"^ over all spatial 
positions at which emission is detected. Detection limits are derived on a cloud-by-cloud basis from 
the rms scatter in emission-free regions of the column density maps. In practice this means that 
cloud boundaries are defined at different absolute emission levels for each object. While less than 
ideal, this data-driven approach is unavoidable given the different noise and intensity levels in the 
three regions. 



-7- 



In mass estimation we assume a distance of 2.6 kpc for clouds A and B and 6.5 kpc for cloud 
C, and include a factor of 1.4 to account for the presence of helium and heavier elements. Masses, 
peak column densities, and detection limits are summarized in table [TJ along with cloud sizes and 
peak main beam temperatures. 



3.2. Cloud A: Active Star Formation &; Translucent Molecular Gas Embedded in 

the Shell Walls 

3.2.1. Morphology and Physical Properties 

Cloud A is the brightest and most massive of the three objects, with estimated masses of 
4.9 x 10 3 M Q , 2.0 x 10 3 M and 460 M Q in the regions traced by the 12 CO, 13 CO and Cho- 
lines, respectively. Figure [4] shows velocity channel maps of all three tracers, revealing much newly- 
resolved structure that was not visible in earlier low-resolution maps. 

The cloud is dominated by a bright nucleus that stretches from around 4.5° to 4.6° in Galactic 
latitude. The material in this region is traced in both 12 CO(J=1-0) and 13 CO(J=1-0), indi- 
cating relatively dense and high column density gas. The two 13 CO emission peaks at (l,b) = 
(285.91°, 4.52°) and (285.90°, 4.57°) each host sizable C 18 cores, with linear extents of ~ 2-3 pc 
and estimated molecular masses of 350 Mq and 110 M & respectively. Peak column densities in both 
cores are Nk 2 ~ 8 x 10 21 cm -2 , and number densities are estimated to be n ~ 1500. The presence 
of so much C 18 implies dense, possibly star-forming molecular gas, and indeed this entire central 
zone is associated with multiple indicators of active star formation, as will be discussed below. 

Above this dense nucleus, from latitudes of b ~ 4.63° to the limits of the map at b ~ 4.88°, lie 
several fainter filamentary features. These extended, wispy structures are only partially detected 
in 13 CO(J=1-0), indicating relatively diffuse, lower column density material. With angular widths 
of ~ l'—4f (~ 1-3 pc), and orientations parallel to the main shell wall, these features are strongly 



reminiscent of substructure observed at higher latitudes in this section of wall in Paper I There we 
noted the presence of thin, twisted filaments between 4.9° < b < 6.2°, with widths as narrow as the 
~ 2.5' resolution limit of the data. These filaments were detected in Hi and also traced partially 
in CO, highlighting the tight relationship between the molecular and atomic material. The present 
data demonstrate that similar structural characteristics persist well into the most molecule-rich 
regions of the wall. In a scenario where Cloud A has formed in-situ from the swept-up medium, 
these diffuse molecular filaments are convincingly interpreted as material that is less advanced along 
an evolutionary pathway that runs from atomic to molecular to star- forming gas. 
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3.2.2. Star Formation 



The bright core of Cloud A is actively forming high mass stars. It is home to the 2MASS 
stellar cluster DBSB 49 ( |Dutra et al.|2003[ ) - a small open cluster with an estimated age of 2.1 ±0.3 



Myr, whose most massive member is of spectral type BOV (Soares et al. 2008). Two bright IRAS 



embedded young stellar object (YSO) candidates, IRAS 10489-5403 and IRAS 10492-5403, also lie 
within this core region, providing further evidence that star formation is still ongoing. Here, the 
IRAS source selection criteria are a 'high' or 'moderate' data quality flag at 25 ^m and 60 /im, 
and a 60/25 fj,m flux density ratio of > 1.0 (see Beichman et al.|[l~9 88). Assuming a spectral energy 
distribution that peaks at 100 /Lim ( Myers et al.||i~9 87), the bolometric luminosities of these sources 
are estimated to be 1200 Lq and 1700 Lq - equivalent to ZAMS spectral types of ~B4-B5 if each 
source corresponds to a single object. The brighter of the two, IRAS 10492-5403, resides in the 
highest column density region of the lower C 18 core 

The cluster is powering an optical Hn region ~ 15' (~ 11 pc) in diameter. This Hn region 
is seen in SHASSA Ha emission maps dGaustad et~al]|2001[ ) as a bright compact source over an 



order of magnitude brighter than the weakly ionized inner rim of the shell walls (see Dawson et al 
2008bp . Figure [5] shows the red DSS2 optical image of this region, overlaid with the positions of 
the stellar cluster and IRAS sources described above, with 13 CO and Ha contours shown on an 
adjoining panel. The lowest Ha contour marks the estimated boundary of the Hn region, and is 
set at the level at which SHASSA Ha emission drops to approximately 1.2 times that of the local 
background. This boundary extends far beyond the edges of the CO emission, indicating that 
the Hn region is no longer confined by its parent molecular cloud. Prominent dark dust lanes at 
(l,b) ~ (285.91°, 4.52°) coincide with the highest column-density regions of the molecular gas, and 
the small molecular clump at (285.87°, 4.39°) is also seen in absorption against the optical nebula. 

The ionizing flux required to power the Hn region provides an independent probe of the 
stellar population of the cluster. The expression for the radius of a classical Stromgren sphere, 
R = (3N* /Airan 2 ,) 1 / 3 , may be used to obtain a first order estimate of this flux under the assumption 
that ionization equilibrium has been reached. We estimate the electron density, n e , from the 
emission measure via the relation, EM (cm -6 pc) = 2.75 T®' 9 In a e 2 ' 2E<yB ~ v " > = n 2 L, where T4 
is the temperature in 10 4 K, Iu a is the Ha intensity in Rayleighs, E(B — V) is the reddening 
due to material in front of the Hn region, and L is the linear depth of the emitting column. The 
reddening is obtained from the neutral gas column density via the relation E(B — V) = ^4y/3.1 = 
A^h/(5.9 x 10 21 cm -2 ) (e.g. Draine & Bertoldi 1996). For the peak Ha intensity of 140 R, a linear 

where the lower 



size of L = 11 pc, and a typical temperature of 10 K, we obtain n e ~ 6-15 cm" 
value arises from ignoring the reddening term completely, and the higher from assuming the full 
column of molecular material at the peak position is located in front of the ionized gas. Taking 
the radius of the Hn region as 5.5 pc, N* = [0.22-1.37] x 10 48 s . This value can be comfortably 



provided by the brightest cluster member, the BOV star (Vacca et al. 1996). We therefore find no 
need for additional massive stars hidden within the cloud. 
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Nevertheless, the level of star formation in Cloud A is highly unusual for such a small, high 
altitude molecular cloud. The scale height of massive star forming regions in the Galactic Disk 
is ~ 30 pc, very similar to the scale height of main sequence OB stars (Bronfman et al. 2000 



Urquhart et al. 


2011 


Snell et al. 


2002), at , 



200 pc the star forming nucleus of Cloud A is still an extreme outlier in 
this distribution. In addition, the cloud mass is somewhat lower than typically observed for objects 
hosting B0 stars (Larson 1982 Dobashi et al. 2001), suggesting a high star formation efficiency. 
The unambiguous identification of Cloud A as part of an expanding superstructure highlights the 
likely importance of supershells in providing star-forming molecular gas to the upper regions of the 
Galactic Disk. 



3.2.3. Dynamics 

The star formation activity strongly influences the dynamics of Cloud A and the region sur- 
rounding it. The molecular gas in the bright nucleus shows wide or multi-component spectral 
profiles, with typical velocity dispersions of Av ~ 2.35a v ~ 2-3 km s _1 in 12 CO(J=1-0); signifi- 
cantly broader than elsewhere in the cloud. This dynamical influence is even more pronounced in 
the structure of the atomic material in the shell wall. Figure [6] shows Hi channel maps overlaid with 
12 CO (J= 1-0) and SHASSA Ha contours. The Hn region has carved out a cavity in the neutral 
gas, about which the Hi shows the classic signature of expansion - a wide void at vi sr ~ —23 km 
s _1 that fills in at the velocity channels to either side. This behavior manifests in the Hi spectra as 
a characteristic broadened or double-peaked profiles, some examples of which are shown in Figure 
[7j This cavity opens out to the shell interior at its lower extreme, indicating that the Hn region 
has not only broken out of its parent molecular cloud, but is also in the process of blistering out of 
the atomic shell wall in which the whole complex resides. 

Because the presence in Hi spectra of two peaks flanking a trough is also a characteristic 
signature of absorption - either against a continuum background source or Hi self absorption 
(HISA) - we take a moment to examine these two alternative hypotheses. In the case of continuum 
absorption, the 1.4 GHz continuum brightness of the Hll region, as derived from the continuum 
data present in the original Hi observations, is only < 2 K - too low by an order of magnitude 
to explain the observed depths of the troughs. In the case of HISA, even when a clear trough is 
present, the outer slopes of the lines are too steep to be fitted by a single positive warm background 
component and negative cold foreground component, and are instead well fitted by two positive 



Gaussians (see also Paper I). Most tellingly, distinct double-peaked or multi-component spectra 
are observed in all three CO lines over the same region. Unlike Hi these profile shapes cannot 
arise as a result of absorption, and must indicate the presence of separate velocity components. 
Taken together with the spatio-velocity structure of the neutral gas described above, the evidence 
therefore strongly favors a dynamical origin for the split profiles observed throughout this region. 



Despite a strong correlation between the spatial positions of multi-component Hi and CO 
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spectra, the relationship between the emission components in the two tracers is generally not 
simple. Although both Hi and CO may exhibit multiple velocity components at any given position, 
these components do not in general correspond exactly, indicating some independence of movement 
between the atomic and molecular material, and more complex dynamics than a simple expanding 
shell. In particular, the velocity separation between CO components is typically no greater than 
~ 2.5 km s , whereas for Hi this value reaches ~ 5.5 km s -1 , with the largest separations observed 
in regions where CO is not detected. The origin of this behavior may lie in the different densities of 
the molecular and atomic components, the former of which is less easily accelerated by the energy 
input from the star forming region. 

The energy input that drives the expanding motion of the ISM around massive stars is provided 
both by over pressurization of the warm Hn region, and by the mechanical energy input of stellar 
winds. Combined Hn region and stellar wind models predict that ~ 0.1% of the combined energy 
input into the system in the form of H-ionizing photons and wind luminosity is converted to kinetic 



energy of the surrounding gas over the lifetime of the star (Freyer et al. 2003, 2006). Assuming 



separate expansion velocities of 2.75 km s 1 and 1.25 km s 1 for the atomic and molecular compo- 
nents, and associated masses of ~ 900 M and ~ 3.0 x 10 3 M , respectively, the kinetic energy of 



the neutral gas is estimated to be ~ 1.1 x 10 47 ergs. Here, the atomic mass is estimated from the 



fitted Hi spectra of Paper I, assuming an optically thin medium, and the molecular mass is derived 
from 12 CO(J=1-0) emission, with both tracers summed over the vicinity of the Ha nebula. A B0V 
star emits Lyman continuum radiation at a rate of ~ 1.4 x 10 48 photons s _1 (Vacca et al. 1996), 
which gives a lower limit of ~ 2.0 x 10 51 ergs for the total energy input in the form of ionizing 
radiation over the 2 Myr lifetime of the cluster, assuming a wavelength of 91.2 nm for all photons. 
The contribution from the stellar wind is considerably lower. The mass loss rates and terminal 
velocities of Schaerer Sz de Koter (1997) imply a total mechanical energy input of only ~ 7.6 x 10 48 



ergs for a B0V star over the cluster lifetime. The total energy input in the form of Lyman con- 
tinuum photons and mechanical wind luminosity is therefore taken to be ~ 2.0 x 10 51 ergs. The 
measured kinetic energy of the neutral gas is only ~ 0.01% of this figure, implying that the energy 
input of the star is at least an order of magnitude greater than that required to drive the motion of 
the surrounding ISM. This discrepancy may arise partially from an underestimate of the associated 
neutral gas mass. Over half of the material in this section of wall is known to be either cold Hi or 
CO-dark H2, and is hence not recovered by the standard mass estimation methods employed here 



(see Paper I ) . Secondly, the fraction of input energy lost to the system by dust heating and cooling 



may be significant, and is not accounted for in the models of Freyer et al. (2003, 2006). Moreover, 



unlike the idealized model case, the Hn region has broken out of its confining section of shell wall, 
allowing the escape of energy into the shell interior. Taking these factors into consideration, we 
consider the energy input of the central star to provide a reasonable explanation for the dynamics 
of the gas of the star-forming region. 

Outside the influence of the star forming region, above about b ~ 4.65°, CO velocity dispersions 
are smaller (Av ~ 2.35a v ~ 1-2 km s _1 in the 12 CO line) and spectral profiles simpler. In cases 
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where double components are present, they can generally be assigned to discrete spatial features in a 
way that is not possible in the more confused emission of the star-forming zone. Hi spectra are also 
generally well fit by single Gaussian profiles, although the presence of some barely resolved double- 
component lines (see Paper I[ ) hints at further complexity that might be more clearly revealed by 
data at higher spatial or spectral resolutions. Nevertheless, there is no evidence to suggest the kind 
of deterministic motions observed in the lower regions of the cloud, and velocity structure here is 
well explained as a combination of turbulent and bulk motions imprinted on the wall during its 
formation. This difference between the dynamics of the star-forming nucleus and the more diffuse 
material outside it highlights a transition between a zone in which the properties of the ISM are 
governed by the physical conditions in the superstructure of which it is part, and the stage at 
which the local influence of star formation starts to dominate. Current estimates of the ages of 
GSH 287+0417 (~ 10 7 yr) and the cluster (2.1 ± 0.3 x 10 6 yr) imply that this transition began 
relatively recently in the shell's lifetime. 



3.3. Clouds B & C: Dynamical Disruption of Pre-existing Molecular Gas & 

Long-Lived Remnant Clumps 

3.3.1. Morphology and Physical Properties 

Cloud B is smaller and less massive than Cloud A, with estimated masses of 830 Mq, 310 Mq 
and 80 M traced in 12 CO(J=1-0), 13 CO(J=b-0) and C 18 O(J=l-0), respectively. However, the 
fraction of the total mass traced by the rarer 13 CO and C 18 lines is the same for both objects, at 
approximately 40% and 10%. 

Cloud C is the smallest of the three clouds, with estimated masses 350 Mq and 70 Mq traced 
in 12 CO(J=1-0) and 13 CO(J=1-0), respectively. The mass fraction traced by the 13 CO line is thus 
approximately half that of the other two clouds. No C 18 O(J=l-0) is detected, to an estimated 
H2 column density detection limit of 1.2 x 10 21 cm -2 . This translates to a mass detection limit of 
~ 40 Mq, meaning that a C 18 core with a fractional mass of ~ 10%, as seen in Clouds A and B, 
would not be detected in the present data. 

Figures § and § show 12 CO(J=1-0), 13 CO(J=1-0) and C 18 O(J=l-0) velocity channel maps of 
Cloud B, and the relationship of this molecular gas to the surrounding and atomic material. The 
data reveal new morphological evidence to support the interpretation that this dense-tipped head- 
tail feature was formed by the Between (I, b) ps (288.30°, 5.58°) and (288.24°, 5.69°) the cloud shows 
a characteristic bow-like shape, with a sharp intensity gradient along the curved front edge facing 
the centre of the shell. This sharp edge is seen in both 12 CO(J=1-0) and 13 CO(J=1-0), with the 
latter concentrated towards the front of the feature, suggesting the preferential build up of denser 
gas there (see also Figure [3]). This bow-shaped zone also hosts a C 18 core with an estimated mass 
of ~ 80 Mq and density of n ~ 1500 cm -3 , located at the position of peak 13 CO(J=1-0) intensity. 
To the right, the edge of the molecular cloud neatly hugs the border of the Hi emission, tapering 
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to a narrow point at (l,b) = (288.24°, 5.69°). Behind this front, trailing fragments are seen at at 
(I, b) ~ (288.31°, 5.65°), following the direction of the Hi tail. These features are predominantly 
detected in 12 CO alone, and appear to represent more diffuse molecular gas breaking away from 
the back edge of the cloud. 

In contrast, a bright 13 CO-rich clump is seen lagging behind the expansion at (I, b) ~ (288.31°, 5.50°) 
At the present resolution of the Hi data, the peak of this molecular clump lies on the very edge 
of the atomic gas feature and may represent a particularly dense fragment of the disrupted cloud 
left behind as the surrounding envelope is stripped. The emission joining this clump to the rest of 
the cloud shows an unusually straight edge, and very square corners where it intersects the rest of 
the cloud. Surprisingly, this very regular squarish shape accurately reflects the global morphology 
of the gas. This can be readily verified by reprocessing the data from different scan directions 
separately, which demonstrates that the 90° angle between the main components of the cloud is 
real. However, the extreme squareness is partially enhanced by the coincidental alignment of cloud 
structure with the telescope scan direction. 

Cloud C, at a relatively large distance of 6.5 kpc, is resolved here for the first time, and shows 
many morphological similarities to its nearer cousin. Figures M and 10 show 12 CO (J= 1-0) and 



13 CO(J=1-0) velocity channel maps, and the relationship between this molecular material and the 
nearby atomic shell wall. The CO emission shows a distinct head-tail structure, with a bright peak 
located in both tracers at (I, b) » (278.42°, 1.00°), and weaker emission extended to the upper-left, 
pointing away from the centre of the shell. The head lies at the very tip of an atomic protrusion, 
about which the curved walls of the shell appear to have been carved out. As in the case of Cloud 
B, the edge facing the centre of the shell shows a sharper intensity gradient in both 12 CO and 
13 CO than elsewhere in the cloud, with 13 CO emission concentrated towards the leading edge. The 
molecular tail shows a second weaker emission peak at (l,b) ~ (278.44°, 1.02°), and beyond that a 
single drop of 12 CO, trailing behind the main cloud at (I, b) « (278.46°, 1.05°). This barely detected 
feature is reminiscent of the trailing fragments observed in Cloud B. 

The wall of GSH 277+00+36 shows several examples of similar Hi drips, spaced periodically 
along the shell walls, with molecular matter so far detected at the tips of only two. Their ori- 
gin has so far not been conclusively determined, and existing works have advanced the differing 
(although not necessarily mutually exclusive) interpretations of these features as Rayleigh Taylor 
drips dMcClure-Griffiths et aJ.|2003| | or as protrusions left behind through interaction of shell with 



pre-existing dense clumps ( Paper I ) . The present observations provide further support for this latter 
interpretation. 



3.3.2. Cloud Dynamics and Clump Survival 



The spectral profiles of both CO and Hi in Clouds B and C are simple, and generally well 
fit by single Gaussian components at coincident velocities. 12 CO(J=1-0) linewidths are typically 
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Avfwhm ~ 1.0—1.5 km s -1 , and we see no direct kinematical evidence of cloud-shell interaction, 
such as shock-broadened line wings. This latter point is unsurprising. Based on the expansion 
velocities and projected distances of the molecular clouds from the main walls, the shell-cloud 
interaction would have occurred > 4 Myr in the past for Cloud B and > 1 Myr in the past for 
Cloud C. Assuming the velocity of the shock in the cloud material is given by v c i = v s h (p s h/ Pel) ' 5 
(e.g. Klein et al"|l994 ), where v s h is the shell expansion velocity and p s h/p c i is the mean shell-cloud 
density contrast, shocks will propagate completely through Clouds B and C on similar timescales 
of ~ 4 Myr and ~ 1 Myr, respectively, for shell wall densities of = 10 cm -3 (Paper I) and 
mean cloud densities of nu = 600 cm~ 3 and 250 cm -3 , as estimated from present 12 CO data. 
This suggests that the shell shock fronts have now passed through both clouds, and initial shock 
signatures will have begun to dissipate. Any current dynamical interaction with the wider shell 
system is now presumably limited to the flow of the low-density interior medium past the remaining 
molecular gas; an interaction that is unlikely to produce strong spectral signatures. Moreover, 
shocked molecular gas is generally more readily traced in the higher excitation CO rotational lines, 
and may remain invisible in CO(J=1-0), even when present (e.g. Reach et al. 2005). 

As pointed out in Paper I[ extinction arising from material outside the CO gas in Cloud B is 
insufficient to shield it against FUV (6 eV < hv < 13.6 eV) photo-dissociation, and the transition 
between the molecular and atomic phases of the ISM is presumably proceeding in the molecular 
to atomic direction. However, while CO close to the surface of the cloud will be quickly destroyed, 
the high column densities at its peak positions imply long survival times for molecules deep within 
it. The CO destruction rate may be expressed as 

^CO^o/co e~ 2bcoAv s" 1 (see Appendix C of 
Wolfire et al. 2010), where I' co is the unshielded CO photodissociation rate in the Draine field, 



G' is the ratio of the incident radiation field to the Draine field, fco is the CO self shielding 
factor, and bco is the dust extinction coefficient for CO extinction. Estimating Ay from Nr 2 
using the relation Ay = Nn/(1.9 x 10 21 cm -2 ), and taking all other parameters as in Paper I, we 
obtain characteristic dissociation timescales of r > 10 s yr for molecules at the centre of the bright 
peaks of the cloud. While this simplistic treatment does not account for the time-variation of the 
destruction rate with decreased self-shielding nor the ongoing recombination of CO molecules, it 
nevertheless suggests that the gas in its present state cannot be easily rendered unobservable by 
photodissociation, provided the background UV field is not substantially higher than in the solar 
neighborhood. This is a reasonable assumption given the low level of ionization observed around 
the inner rim of GSH 287+04-17 pawson et al"1|2008b[ ). 

Other destructive influences on the molecular gas include thermal evaporation from the hot in- 
terior medium and dynamical interaction with the shell system. The classical evaporation timescale 
for Cloud B was estimated to be ~ 10 s yr (Paper I), and this figure does not change substantially as 
a result of the resolution improvements of the present observations. In contrast, dynamical interac- 
tion with the shell is likely to initially disrupt the clouds on much shorter characteristic timescales 
of a few Myr (Paper I). However, the most intense phase of this interaction is limited to the finite 
period in which the initial shock imparted through interaction with the dense material of the shell 
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is still propagating through the cloud; a stage that Cloud B has already survived. For a canonical 
interior medium density of ~ 0.01 cm -3 and the slow flow-speed associated with the expansion of 
the shell, the current dynamical influence on the clouds is minimal. 

In the absence of a continuing disruptive influence from the shell system, the fate of Cloud B 
will be determined by the balance between its own internal velocity dispersion and the confining 
forces of external pressure and gravity. Assuming a spherical cloud with a r _1 radial density profile 
and ignoring external pressure and magnetic fields, the virial mass required for self-gravity to just 
balance the system is given by M v i r /M & = 190 R (Av c id) 2 , where R is in parsecs and Avdd is the 
FWHM velocity dispersion in km s _1 . Under this formulation, although the cloud as a whole is 
marginally unbound, the virial mass of the C 18 core is just sufficient for it to be self-gravitating 
(M vir = 80 M Q ). 

Performing the same analysis for Cloud C, we may make an identical argument concerning the 
minimal destructive influence due to evaporation by the internal medium and the current dynam- 
ical interaction of the cloud with the shell system. However, the characteristic photodissociation 
timescale for CO at the centre of the cloud is found to be much shorter, at ~ 10 6 yr. This timescale 
is very sensitive to beam dilution effects, however, and may increase significantly in better resolved 
data. In addition, Ha data for this region is confused, precluding comment on the likely magnitude 
of background UV field, and rendering the analysis inconclusive. The virial mass of Cloud C is 
found to be several times larger than its CO intensity derived mass, suggesting that it will disperse 
if not confined by external pressure. 



3.3.3. Comments on Star Formation 



Self-gravitating C 18 cores are a good predictor of star-forming gas (Tachihara et al. 2002), 
and we might therefore hope to discover signs of ongoing star formation in Cloud B. However, 
although several IRAS sources fall within the cloud boundary, no clear evidence is found. A single 
IRAS source with YSO colors is located at (l,b) = (288.29°, 5.58°), but is not associated with any 
other known star formation tracers or major CO intensity peaks. Two sources that do not fulfill the 
YSO selection criteria are also located within the cloud boundaries at (I, b) = (288.26°, 5.59°) and 
(288.30°, 5.50°). These sources are detected at 100 //m and 60 fim and at 60 fim only, and are both 
located close to CO emission peaks. However, in the absence of any corroborating evidence, it is 
unclear whether any of these sources trace embedded star formation or point-source-scale structure 
in the cloud itself. Nevertheless, the low bolometric luminosities of all three (L& D z < 100 Lq) allow 
us to constrain any star formation that may be occurring to low mass objects. 
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4. Summary & Conclusions 

The transition between atomic and molecular gas forms a key link in the lifecycle of material 
in galactic systems. In a dynamic and multi-phase ISM, large-scale stellar feedback can drive this 
transition in either direction. Diffuse material accumulated and compressed in supershells can reach 
sufficient densities and column densities for the formation of molecular gas. Conversely, pre-existing 
molecular clouds may be disrupted by their encounter with an expanding shell shock front, leaving 
them stripped of shielding material and vulnerable to UV dissociation. 

In this context we have performed an in-depth case-study of three molecular clouds in the 
walls of the Galactic supershells GSH 287+04-17 and GSH 277+00+36. These clouds have been 
identified in previous work as likely locations at which the atomic-molecular transition is proceeding 
in opposite directions. All three objects were mapped to a resolution of ~ 45" in the 12 CO(J=l- 
0), 13 CO(J=1-0) and C 18 O(J=l-0) lines, providing detailed information on the distribution and 
properties of the molecular ISM, and enabling a greatly improved discussion of its relationship to 
the wider environment in which it resides. 

Clouds B and C are small condensations of molecular gas with molecular masses of ~ 830 M & 
and ~ 350 M , respectively, whose location at the tips of finger-like extensions of Hi suggests the 
sweep-over and subsequent entraining of pre-existing clouds. Our observations reveal hitherto un- 
seen morphological signatures of disruption and possible shock-interaction in both clouds, strongly 
supporting this interpretation. Both objects show dense heads with sharp intensity gradients at 
their leading edges, and trailing tails of more diffuse material. The emission of Cloud B, the better 
resolved of the two, closely delineates the curved front rim of the long Hi feature with which it is 
associated, with molecular gas forming a bow-like shape - a characteristic signature of a shocked 
cloud. 

Both clouds have survived the original interaction with their respective shells, and are now 
interacting with the diffuse interior medium. High density contrasts between the molecular gas 
and this interior medium suggest that the destructive dynamical and thermal effects of the shell 
systems at the present epoch are minimal; i.e. while the interaction with an expanding shell may 
have carved the ISM into the configurations we currently observe, the shell systems are no longer 
a primary determinant of the cloud lifetimes. In the case of Cloud B, high peak column densities 
and an absence of nearby sources of FUV photons implies that although photodissociation of CO 
is underway, and the gas is becoming atomic, this process is proceeding at a relatively gentle 
pace, suggesting potentially long survival times if the cloud remains in its current configuration. 
However, the internal velocity dispersions of both Clouds B and C suggest that neither are globally 
self-gravitating, and may disperse if not confined by external pressure. Cloud B does, however, 
contain a dense C 18 core which is locally bound. Such cores are generally good predictors of 
star-forming gas, and the presence of several cool, low-luminosity IRAS sources within the cloud 
hints at possible low mass star formation activity. 

Cloud A is a moderately sized, massive-star-forming cloud embedded within the main wall of 
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GSH 287+04-17. With a substantial molecular mass of ~ 4.9 x 10 3 M , and a distribution that 
closely follows that of the curved atomic shell, it has been identified in previous work as a likely 
example of molecular material formed in-situ from the swept-up atomic ISM. Our observations 
reveal a bright star-forming nucleus detected in all three CO lines, from which emanate several 
diffuse filaments detected primarily in 12 CO(J=1-0) alone. 

The star formation in Cloud A is a convincing example of second generation massive star 
formation triggered by the primary generation that formed the shell. The star forming nucleus 
hosts the young (2.1 ± 0.3 Myr) cluster DSBS 49, whose most massive member is of spectral type 
BOV. The formation of such massive stars is highly unusual at this Galactic altitude (z ~ 200 pc), 
and the unambiguous identification of this region as part of an expanding superstructure highlights 
the role of supershells in providing star-forming material to the upper regions of the Galactic Disk. 
The cluster is powering an optical Hn region ~ 11 pc in diameter, which has carved out a cavity in 
the shell wall and is in the process of blistering out into the interior. Expansive motions driven by 
this star formation activity dominate the dynamics of both the atomic and molecular gas in this 
local region. 

In contrast, the filaments outside of the influence of this star- forming zone form quiescent parts 
of the shell wall. These diffuse, wispy-looking molecular features are ~ 1-3 pc in width, and closely 
resemble twisted, filamentary structure detected in both Hi and CO elsewhere in the same section of 
wall (Paper 1), demonstrating that common structural characteristics persist throughout much of 
its volume. In a scenario where Cloud A has formed in-situ from the swept-up medium, these diffuse 
molecular filaments are interpreted as material that is less advanced along an evolutionary pathway 
that runs from atomic to molecular to star-forming gas. The border between these filaments and 
the star forming nucleus of the cloud marks a transition in which the properties of the ISM switch 
from being determined by the global characteristics of the swept-up material to being dominated 
by the local influence of the stellar sources forming within it. This critical stage in the process of 
triggered star/cloud formation is currently ongoing in the walls of GSH 287+04-17, ~ 10 yr after 
its birth. 

The data presented in this work mark all three clouds as excellent targets for further obser- 
vation, setting the stage for an empirical study of the transition both ways between molecular and 
atomic gas. 

The Australia Telescope Compact Array, Parkes and Mopra Telescopes are part of the Aus- 
tralia Telescope which is funded by the Commonwealth of Australia for operation as a National 
Facility managed by CSIRO. We acknowledge the use of the Southern H- Alpha Sky Survey Atlas 
(SHASSA), which is supported by the National Science Foundation. The DSS2 image used in this 
work is based on photographic data obtained using The UK Schmidt Telescope. The UK Schmidt 
Telescope was operated by the Royal Observatory Edinburgh, with funding from the UK Sci- 
ence and Engineering Research Council, until 1988 June, and thereafter by the Anglo-Australian 
Observatory. Original plate material is copyright of the Royal Observatory Edinburgh and the 
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Anglo- Australian Observatory. The plates were processed into the present compressed digital form 
with their permission. The Digitized Sky Survey was produced at the Space Telescope Science 
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Table 1. Summary of Cloud Properties 



12 CO(J=1-0) 13 CO(J=1-0) C 18 O(J=l-0) 



Cloud T max a A r H 2 , ii m h WH 2 ,mni C ^ci/ Mdd c T ma x^ ^h 2 ,Um h Nn 2 ^ max c j4 c ; £ j d M c ;,j e T max a Nn 2 j im b Nu 2imax c M c id" 
K 10 21 cm- 2 10 21 cm- 2 pc 2 Mq K 10 21 cm" 2 10 21 cm- 2 pc 2 M K 10 21 cm- 2 10 21 c m - 2 pc 2 M Q 



A 21.8 0.56 9.0 155 4900 7.4 0.55 9.9 57 1950 1.1 3.6 8.5 4.8 460 

B 12.5 0.32 2.7 43 830 4.3 0.40 2.2 18 310 1.1 2.4 5.1 1.1 80 

C 5.4 0.24 1.1 39 350 1.0 0.15 0.46 18 70 < 0.2 1.2 ■■■ ■ ■ ■ < 40 

a Peak main beam temperature 

b 3ir detection limit in H2 column density. 

c Peak value of H2 column density. 

d Combined area of all pixels above H2 column density detection limit. 

c Total molecular gas mass traced in each CO line. Includes factor of 1.4 to account for presence of He and heavier elements. 
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Fig. 1. — Integrated intensity image of a subsection of the wall of GSH 287+04-17. The greyscale 
image is Hi observed with the ATCA and Parkes telescopes, and contours are 12 CO(J=1-0) observed 
with the NANTEN telescope (both lifted from Dawson et al.|[20TT ) . The velocity integration range 
is —26.5 < vi sr < —19.9 km s , and CO contour levels are drawn at 1.5+5.0 K km s _1 . Dashed 
lines show regions mapped with the Mopra telescope in the present work. 
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Fig. 2. — Integrated intensity image of a subsection of the wall of GSH 277+00+36. The greyscale 
image is Hi observed with the ATCA and Parkes telescopes, and contours are 12 CO(J=1-0) observed 
with the NANTEN telescope (both lifted from Dawson et al.|[20TT ) . The velocity integration range 
is 41.6 < vi ST < 44.1 km s _1 , and CO contour levels are drawn at 1.2+1.0 K km s _1 . Dashed lines 
show regions mapped with the Mopra telescope in the present work. 



- 21 - 
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Fig. 3. — Nu 2 maps of Clouds A, B and C. Color images show material traced in 12 CO(J=1-0), 
black contours in 13 CO(J=1-0) and white contours in C 18 O(J=l-0). Contour levels begin at the 3a 
detection limit for each tracer and are as follows. Cloud A - 13 CO: [0.55, 1.05, 1.55, 2.55, 3.55, 4.55, 
5.55, 6.55, 7.55, 8.55] xlO 21 cm" 2 ; C 18 0: [4.5, 5.5, 6.5, 7.5] xlO 21 cm- 2 . Cloud B - 13 CO: [0.40, 
0.65, 0.9, 1.15, 1.40, 1.65, 1.9]xl0 21 cm" 2 ; C 18 0: [2.5, 3.5, 4.5]xl0 21 cnr 2 . Cloud C - 13 CO: [0.15, 
0.25, 0.35, 0.45] xlO 21 cm" 2 . In cloud A, C 18 contours are only plotted within the 12 CO(J=1-0) 
cloud boundary (i.e. noise peaks from high- noise regions outside the cloud have been removed). 
Crosses show IRAS point sources with moderate or high quality flux measurements at 60 /jm and 
25 ^m, and Fqq > F25. 
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Fig. 4. — Velocity channel maps of the embedded cloud, Cloud A. Each panel shows an average 
over four channels, corresponding to an interval of 0.7 km s _1 . The color scale shows 12 CO(J=1-0) 
emission, black contours show 13 CO(J=1-0) and white contours show C 18 O(J=l-0). Contour levels 
start at the 4<r detection limit for each dataset, corresponding to 0.60 and 0.45 K for 13 CO and 
C 18 respectively. Contour intervals in 13 CO are every 0.5 K until 1.6K and every 1.0 K thereafter. 
Contour intervals in C 18 are every 0.1 K. C 18 contours are only plotted within the 12 CO(J=1-0) 
cloud boundary (i.e. noise peaks from high- noise regions outside the cloud have been removed). 
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Fig. 5. — Anatomy of the star forming region in Cloud A. The color image shows optical data 
from the second generation Digital Sky Survey (DSS2). The dashed circle marks the radius of 
the young IR cluster DBSB 49, while red crosses mark IRAS point sources with YSO colors. 
The contour image shows velocity-integrated 13 CO(J=1-0) emission (black contours), scaled to 
units of H2 column density, with levels at [0.7, 1.7, 2.7, 4.7, 6.7, 8.7] xlO 21 cm -2 . Pink contours 
are SHASSA Ha emission beginning at a level of 55 Rayleigh (R) and incremented every 20 R 
thereafter. Numbered blue circles mark the locations of the spectra shown in Figure [7j 
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Fig. 6. — Velocity channel maps of Hi and 12 CO(J=1-0) in the vicinity of Cloud A. Color panels 
are Hi. Black contours are 12 CO(J=1-0), beginning at 2.5 K and incremented every 3.0 K. White 
contours show SHASSA Ha emission (with no velocity information) beginning at a level of 55 
Rayleigh (R) and incremented every 20 R thereafter. The black dashed line marks the limits of the 
region observed in CO. 
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Fig. 7. — Hi and CO(J=1-0) spectra at the positions numbered in Figure [5j Hi data are from 



Dawson et al. (2011) and CO data have been smoothed to match the 2.5' resolution of the Hi 
cube. The left- and right-hand scales show Hi and CO intensities, respectively. 13 CO(J=1-0) and 
C 18 O(J=l-0) spectra have been scaled by factors of 2 and 5 respectively for ease of viewing. 
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Fig. 8. — Velocity channel maps of the offset cloud, Cloud B. Each panel shows an average over 
three channels, corresponding to an interval of 0.5 km s _1 . The color scale shows 12 CO(J=1-0) 
emission, black contours show 13 CO(J=1-0) and white contours show C 18 O(J=l-0). Contour levels 
start at the 4cr detection limit for each dataset, corresponding to 0.7 and 0.4 K for 13 CO and C 18 
respectively. Contour intervals are every 0.5 K in 13 CO and every 0.1 K in C 18 0. 
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Fig. 9. — Hi and 12 CO (J= 1-0) integrated intensity in the vicinity of Clouds B and C. Color is 
Hi and black contours are 12 CO(J=1-0). For the Cloud B region, emission is integrated over the 
velocity range —25.2 < vi sr < —21.9 km s -1 , and CO contours begin at 3.0 K km s~ 1 and are 
incremented every 2.5 K km s _1 . For the Cloud C region, emission is integrated over the velocity 
range 41.2 < V[ sr < 43.7 km s _1 , and CO contours begin at 2.0 K km s _1 and are incremented 
every 1.0 K km s _1 . 




GALACTIC Longitude GALACTIC Longitude GALACTIC Longitude GALACTIC Longitude 



Fig. 10. — Velocity channel maps of the offset cloud, Cloud C. Each panel shows an average over 
three channels, corresponding to an interval of 0.5 km s _1 . The color scale shows 12 CO(J=1-0) 
emission and black contours show 13 CO(J=1-0). Contour levels start at the 4cr detection limit of 
and 0.3 K and are incremented every 0.2 K. The white dashed lines mark the region observed in 
13 CO(J=1-0). 



